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ABSTRACT 

For the first time, we have identified photospheric emission lines in the far-UV spectrum of a white dwarf. They were discovered in 
the Far Ultraviolet Spectroscopic Explorer spectrum of the hot (T^ff « 200 000 K) DO white dwarf KPD 0005+5 106 and they stem 
from extremely highly ionized calcium (CaX AA 1137, 1159 A). Their photospheric origin is confirmed by non-LTE line-formation 
calculations. This is the highest ionisation stage of any element ever observed in a stellar photosphere. Calcium has never been detected 
before in any hot white dwarf or central star of planetary nebula. The calcium abundance determination for KPD0005+5106 (1-10 
times solar) is difficult, because the line strengths are rather sensitive to current uncertainties in the knowledge of effective temperature 
and surface gravity. We discuss the possibility that the calcium abundance is much lower than expected from diffusion/levitation 
equilibrium theory. The same emission lines are exhibited by the [WCE]-type central star NGC 2371. Another CaX line pair (AA 1461, 
1504 A) is probably present in a Hubble Space Telescope spectrum of the PGl 159-type central star NGC 246. 
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1. Introduction 

Observations of hot (pre-) white dwarfs with the Far Ultraviolet 
Spectroscopic Explorer {FUSE) have revealed a large number of 
chemical elements that were never detected before in these ob- 
jects. Their abundances can be used either to probe interior pro- 
cesses in previous stellar evolution phases or to test predictions 
from theories for element diffusion and radiative levitation. 

We have recently identified Neviii lines in the hottest 
(reff ; 150 000 K) non-DA (pre-) white dwarfs, i.e. objects of 
spectral type PGl 159, DO, and [WCE] (Werner etal. 2007). 
The discovery of these lines in the hottest known DO white 
dwarf KPD 0005-1-5 106 was particularly surprising, because this 
proves that its effective temperature must be much higher than 
previously thought (200 000 K instead of 120 000 K). 

KPD 0005 -h5 106 was frequently observed by FUSE as a cal- 
ibration target over its entire lifetime. We have co-added all 
available spectra and obtained datasets with very high S/N ra- 
tio. A careful inspection of spectra taken with different detec- 
tors revealed the presence of two hitherto unidentified emission 
lines. While there is still a large number of unidentified absorp- 
tion lines present in FUSE spectra of hot white dwarfs (WDs), 
the discovery of emission features is unique and was completely 
unexpected. In this Letter we identify them as CaX lines and 
present results of non-LTE modeling in order to confirm their 
photospheric origin and to perform an abundance determination. 

We present observations and line identifications in Sect.|2] 
and describe the modeling in Sect.[3] The results from line- 
profile fits are presented in Sect.H] We conclude with Sect.|5] 



* Based on observations made with the NASA-CNES-CSA Far 
Ultraviolet Spectroscopic Explorer. FUSE is operated for NASA by the 
Johns Hopkins University under NASA contract NAS5-32985. 



2. Observations and Cax line identifications 

The FUSE instrument consists of four independent co-aligned 
telescopes and spectrographs; two with Al-i-LiF optical coatings 
and two with SiC coatings. Taken together, the four channels 
span the wavelength range 904-1 187 A with a typical resolving 
power of /? «i 20 000. Further information on the FUSE mission 
and instrument can be found in Moos et al. (2000) and Sahnow 
etal. (2000). KPD 0005 4-5 106 was observed as a wavelength 
calibration object throughout the FUSE mission (observations 
M1070201-M1070234), in each of the spectrograph apertures 
(LWRS, MDRS, HIRS), and once under program P1040101 in 
LWRS. All observations were obtained in TTAG mode, except 
for three early HIST mode observations; the latter were ex- 
cluded as they have different residual distortions in the wave- 
length scale. All exposures were processed with the final version, 
V 3.2.2, of CalFUSE (Dixon etal. 2007). Spectra for each chan- 
nel were shifted to place absorption lines from low-ionization 
interstellar gas at a heliocentric velocity of -15kms"' (Werner 
etal. 1996), and combined. Variations in losses due to event 
bursts and channel misalignments caused the net exposure times 
to vary: HIRS times varied from 7.0ks in SiC2a to 17.9ks in 
LiF2b; MDRS times varied from 11.6ks in SiCla to 18.6 ks in 
LiF2b; LWRS times varied from 28.5 ks in LiFlb to 33.6 ks in 
LiF2b. The net signal ranges from 700 counts per 0.013 A pixel 
in the HIRS spectra to 1200c/pix in MDRS and as much as 
2700 c/pix in LWRS. The Ca x emission features are seen clearly 
in each of the six available spectra (LiFlb and LiF2a in each of 
HIRS, MDRS, LWRS). These six spectra were then resampled 
onto a common wavelength scale and combined to produce a 
single spectrum with an effective exposure time of 122.2 ks. 

The UV spectrum of NGC 246, taken with the Space 
Telescope Imaging Spectrograph (STIS) with grating E140H 
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Fig. 1. The two CaX emission lines discovered in KPD 0005+5106 (thin graphs). Overplotted is the spectrum from a model with 
Teff - 200 000 K, log g =6.2, and solar Ca abundance. The model was convolved with a Gaussian with FWHM=0.05 A in order to 
match the instrumental resolution. 
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Fig. 2. Grotrian diagram of our Ca x model ion. Lines discussed 
in the text are caused by transitions between n = 4 sublevels. 
The 4p-4d transition causes the observed UV emission lines. 
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Fig. 3. The two Cax emission lines discovered in 
KPD 0005 +5 106 (top spectra) are also seen in the [WCE] 
central star NGC2371. The model profiles are the same 
as in Fig.[T] No attempt is made to fit the possibly wind- 
contaminated [WCE] profiles. For clarity, the FUSE spectra 
were smoothed with Gaussians (FWHM 0.05 and 0.1 A, respec- 
tively). The model spectra (thick lines) were convolved with 
0.1 A Gaussians. 



aboard the Hubble Space Telescope (HST) was retrieved from 
the MAST archive. 

In the FUSE spectrum of KPD 0005+5106 we detected 
two emission lines, located at photospheric rest wavelengths 
AA 1136.5, 1159.2A (Fig.B radial velocity +35kms-i; Werner 
et al. 1996). We identify these fines as due the 4p -4d ^P° tran- 
sition in the Ca X ion (Fig.|2l). Compared to the Ritz wavelengths 
in the NISlQ database, both observed lines are located at wave- 
lengths shorter by 0.3 A. Their NIST ^/-values are \og gi fjk = 
0.23 and 0.46, respectively. The third line component of this 
transition is located at A 116 1.4 A according to NIST, so that 
in reality we expect it to be found at A 1161.1 A. Its ^/-value, 
however, is much smaller (log gifik = -0.47) explaining the fact 
that we cannot detect it in the observation. (This is confirmed by 
our line-formation calculations.) 

We searched for the two CaX lines in other hot DO white 
dwarfs and PG1159 stars, but to no avail. As we will demon- 
strate below (Sect.lU, this is a consequence of the extremely 



http://physics.nist.gov/ PhysRefData/ASD/index.html 



high Teff of KPD 0005 +5 1 06 combined with a relatively low sur- 
face gravity. However, these lines are seen in the very hot, early- 
type Wolf-Rayet central star NGC2371 (Fig.[3]l. For this object 
we did not attempt to fit these lines with our (static) model at- 
mospheres, because the profiles might be affected by the stellar 
wind. 

In the course of our model calculations we found that fur- 
ther CaX lines might be detectable in other wavelength re- 
gions. The 4s^p transition gives rise to a line doublet at 
AA 1461.8, 1503. 8A. Our models predict only marginal emis- 
sion features for KPD 0005+5106, which cannot be detected 
in archival spectra taken with the Faint Object Spectrograph 
aboard HST and high-resolution spectra from the International 
Ultraviolet Explorer. However, for NGC 246 our models predict 
absorption lines that are possibly present in a HST/STIS spec- 
trum (Fig.m. The positions of the tentatively identified absorp- 
tion features in NGC 246 differ from the NIST wavelengths by 
-0.7 and -0.2 A, respectively. 
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Fig. 4. HST/STIS spectrum of the PG1159-type central star 
NGC 246 and computed profiles for the 4s^p doublet of Ca X. 
Their shape reflects the stellar rotation of vsin/ = 70kms 
The observation was smoothed with a Gaussian with 
FWHM=0.03 A. 



From NIST level energies one expects the two strongest 
lines of yet another transition of CaX (4d-4f) to be located 
in the optical UV at AA 3478, 3492 A, respectively. A high- 
resolution spectrum of NGC 246 taken with ESO's Very Large 
Telescope and the UVES spectrograph as part of the SPY sur- 
vey (Napiwotzki et al. 2003) reveals no line features there. Our 
model for NGC 246 predicts absorption lines with a depth of 
only 5% relative to the continuum. The relatively fast rotation 
(vsin; = 70kms"') smears the line features considerably, and 
they remain hidden in the noise. Similar weak absorption pro- 
files for these lines are predicted for KPD 0005+5106, but no 
appropriate observations ai^e available. 



3. Model atmospheres and calcium line-formation 

We have designed a calcium model atom for NLTE line- 
formation calculations. These were performed using and keep- 
ing fixed the physical structure (temperature, densities) of line- 
blanketed NLTE model atmospheres which are described in de- 
tail in Werner etal. (2004). In short, they are plane-parallel and 
in hydrostatic and radiative equilibrium. The models are com- 
posed of He, C, O, and Ne. For KPD 0005+5106, we assumed 
helium-dominated atmospheres with admixtures of C=0.003, 
0=0.0006, Ne=0.01 (mass fractions). The high neon abundance 
was derived from Neviii lines (Werner etal. 2007). The C and 
O abundances are uncertain, because they were derived in ear- 
lier work that assumed that KPD 0005+5106 is relatively cool 
(^eif = 120000K; Werner etal. 1996). We verified that vary- 
ing the C and O abundances within reasonable limits does not 
change the CaX lines significantly. A series of models with vari- 
ous Teff and log g values was computed to study the dependency 
of the Ca X lines on these parameters (see Sect.|4|. For NGC 246 
we adopted Teff = 150000K, log ^ = 5.7, and the composition 
He/C/O/Ne=0.62/0.30/0.06/0.02 (Werner et al. 2007). 

The Ca model atom considers the ionization stages vm-xii, 
represented by 1, 15, 25, 4, 1 NLTE levels, respectively, plus a 
number of LTE levels. In the ions Caix-xi we include 23, 126, 
and 2 line transitions, respectively. Atomic data were taken from 
the NIST, Opacity (Seaton etal. 1994), and IRON (Hummer 
etal. 1993) Projects databases (TIPTOPbas^. Fine-structure 
splitting is accounted for in the final formal solution for the 
synthetic line-profile computation, distributing the level popula- 
tions among sublevels assuming LTE. For all lines we assumed 
quadratic Stark broadening for the profile calculation. 
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Fig. 5. Ionization fraction of calcium as a function of atmo- 
spheric depth in the model with Teflf - 200 000 K, log g - 6.2, 
and solar Ca abundance. 



Particularly for the CaX AA 1137, 1159 A lines the val- 
ues of the oscillator strengths differ between the OP and NIST 
databases. We use the OP values for our NLTE level popula- 
tion iterations, because they are complete, in contrast to the 
NIST database. For the final line-profile calculation we prefer 
the NIST values, which are higher than the OP values by w 25%. 
The diff'erences do not affect our conclusions. 

Photoionization cross-sections are taken from the Opacity 
Project database when available or, otherwise, computed in a 
hydrogen-like approximation. Electron collisional rates were 
calculated with the usual approximation formulae. The Ca model 
atoms that were used for this analysis have been developed in 
the framework of the German Astrophysical Virtual Observatory 
(GAVCQ) project and are provided within the Tubingen Model- 
Atom Database TMAZfl 

4. Results 

Figure [T] shows a fit to the CaX emission lines in 
KPD 0005+5106 with a model Teif = 200 000 K, log g= 6.2, 
and a solar Ca abundance (log Ca = -4.22, mass fraction; 
Asplund etal. 2005). In Fig.|5]we show the ionization structure 
of Ca throughout this model atmosphere. Within the entire line- 
forming region Ca xi is dominant, followed by Ca x. In order to 
achieve the observed emission strength in a model with this tem- 
perature and Ca abundance, the surface gravity must be that low 
(log g -6.2). This is 0.3 dex lower than what is preferred from 
the He ii line spectrum (Werner et al. 2007). We will show, how- 
ever, that the fit to the Ca X lines can be achieved with more than 
one parameter set. 

The occurrence of this line emission can be understood when 
the non-LTE departure coefficients /?, = n^^^^/n^'^ for the pop- 
ulations n, of the involved atomic levels and the line source func- 
tion Si are inspected (Fig.|6]l. The line source function is deter- 
mined by the ratio of the departure coefficients of the lower and 
upper levels (/, j) and can be written as S i/By - [exp{hvij/kT) - 
^/[ibi/bj)exp(hvij/kT) - 1]), where By is the Planck function. 
An overpopulation of the upper level relative to the lower (i.e. 
S i/By > 1) may lead to line emission. In fact, this condition is 
fulfilled in the line-forming region (Fig.|6]l, although both levels 
are underpopulated (i.e. < 1). 



|http://cdsweb.u-strasbg.fr/topbase/| 



http://www.g-vo.org 

http://astro.uni-tuebingen.de/~rauch/TMAD/TMAD.html 
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Fig. 6. Left: Departure coefficients bi for the 4p and 4d levels 
that cause the CaX emission lines. Right: Ratio of 4p^d line 
source function 5; to Planck function By at line centre. The thick 
horizontal lines near log m -Q denote the line formation region. 
Dashed lines correspond to the LTE case (bi - 1 and Si/ By = 1). 
Model parameters as in Fig.|5l 



We have computed a small grid of models with different Tes 
and log g (representing the uncertainties with which these pa- 
rameters are known) and Ca abundances. The results are pre- 
sented in Fig.|2l Generally, high T^ff and low log g is necessary 
in order to bring this line into emission (left and right upper pan- 
els). If Teff decreases (and/or log g increases), the lines first turn 
from emission into weak absorption features and then disappear 
at about 140 000 K. This explains why other DOs do not exhibit 
these lines: they are significantly cooler and have higher gravi- 
ties compai-ed to KPD 0005+5106. 

The sensitivity of these emission lines to the calcium abun- 
dance is complicated and depends on T^ff and log g of the at- 
mosphere (Fig.|7] left and right lower panels). Increasing the 
Ca abundance over the solar value can strongly increase the 
emission height (model Teff = 200 000 K log g -6) or decrease 
it (model Teff = 200 000 K log g = 6.5). For KPD 0005+5 106 we 
achieved a good fit at Teif = 200 000 K, log g =6.2, and solar Ca 
abundance. As mentioned, however, our previous analysis favors 
a gravity higher by 0.3 dex. Increasing the gravity to log g = 6.5 
makes the emission weaker, but this can be compensated by si- 
multaneously increasing the Ca abundance to 3 times the solar 
value and Teff to 220 000 K. A much higher abundance, as pre- 
dicted by diffusion theory, can be excluded. We have calculated 
models with 70 times solar Ca abundance. The emission line 
peak heights hardly change in the Teff = 200 000 K, log g =6.2 
model when Ca is increased from solar to 70 times solar; how- 
ever, detailed inspection of the relative strength of both lines 
shows that it is not in agreement with the observation. In the 
observation as well as in the « solar Ca abundance models the 
11 59 A emission is stronger than the 1 137 A emission, as can be 
expected from the gf -value ratio. In the 70 times solar models 
the emission strength ratio is reversed, in contrast to the obser- 
vation. 

Concerning PG1159 stars, the behaviour of these Ca 
lines is rather similar and therefore not shown in detail 
here. There are seven objects that are hot enough and for 
which FUSE spectroscopy is available. These are the low- 
gravity central stars of planetary nebulae Kl-16, Longmore4, 
RXJ2117.1+3412, NGC246 (Feff = 140000-170000K, 
log g 5.5-6), the higher-gravity objects PG1520+525 and 
PG1144+005 (reff = 150000K, log ^=6.5-7), as well as 
the pecuHar H1504+65 (Feff = 200000K, log ^ = 8). Model 
calculations were performed for all of these objects with 
solar Ca abundance (diffusion is not at work in these objects' 
atmospheres; see Unglaub & Bues 2000). They exhibit the 
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Fig. 7. Profile shapes of the CaX AA 1137, 1159A lines as a 
function of Teft , log ^ , and Ca abundance, as given by the labels. 



Cax AA 1137, 1159 A lines as weak absorption features with a 
maximum depth of 5% of the continuum flux. Such weak lines 
cannot be detected in the available FUSE spectra. 

5. Summary and conclusions 

We have discovered CaX emission lines in far-UV spectra of 
the DO white dwarf KPD 0005+5106. This is the first detec- 
tion of photospheric emission lines in this spectral range of any 
hot (pre-) white dwarf. Provencal et al. (2005) discovered low- 
ionisation emission lines in HST/STIS UV spectra of two rela- 
tively cool (Teff ~ 12 500 K) He-rich white dwarfs (spectral type 
DQ). It was shown, however, that they are chromospheric in ori- 
gin. The CaX lines are the highest ionisation stage of any el- 
ement identified in any stellar photosphere. Our discovery also 
represents the first identification of calcium in hot (pre-) WDs. 

The Ca abundance in KPD 0005+5106 is in the range «;1-10 
times solar A more precise determination from the emission 
lines is not possible. A comparison of this result with predic- 
tions from radiative levitation/gravitational diffusion equilibrium 
theory is difficult because T^ff and log g of KPD 0005+5106 
are outside of the range considered by Chayer etal. (1995; their 
Fig. 20). For the closest parameters (Teif = 130 000 K, log g =7) 
a huge overabundance is predicted (2500 times solar). Our es- 
timate for log g is smaller (6.2-6.5) which would result in an 
even higher overabundance. On the other hand it is impossible to 
make a solid estimate for the effect of the higher Jeff (200 000- 
220 000 K) on the behaviour of the Ca equilibrium abundance, 
because the dominant ionisation stage in KPD 0005 +5 106 is 
Caxi, while it is Caviii in the hottest Chayer etal. model 
(TefF = 130 000 K, log g =7.5). Looking at the behaviour of other 
elements (S, Ar), namely how their equilibrium abundance 
changes when their (respective isoelectronic) ionisation stages 
increase (with increasing Teff ), it is suggestive that the Ca abun- 
dance at Teff =200 000 K is lower than at 130 000 K, but not by 
orders of magnitude. Although detailed calculations are required 
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for a definitive statement, we conclude that the atmosphere of 
KPD 0005+5 106 is probably not in levitation/diffusion equiUb- 
rium. This is confirmed by the diffiision/mass-loss calculations 
of Unglaub & Bues (2000) which suggest that KPD 0005 +5 106 
has yet to cross the wind-limit on its evolutionary track, mean- 
ing that mass-loss is large enough to prevent both gravitational 
settling and the accumulation of radiatively supported heavy el- 
ements. In this case, KPD 0005+5106 is not a descendant of the 
PG1159 stars. An evolutionary link to the He-dominated cen- 
tral stars of spectral type 0(He) and to the RCrB stars has been 
suggested (Werner et al. 2008). 

If unaffected by diffusion processes, then the photospheric 
composition of KPD 0005+5106 is the consequence of previ- 
ous evolutionary phases. In contrast, the presence of Ca in the 
atmospheres of cooler white dwarfs (spectral types DAZ and 
DBZ, with low-ionisation optical Ca absorption lines) requires 
on-going accretion of circumstellar matter, because gravitational 
settling rapidly removes heavy elements from the photosphere 
(e.g. Koester & Wilken 2006). 

The non-detection of the Cax AA 1137, 1159A Unes in the 
hottest PGl 159 stars is explained by undetectably weak absorp- 
tion line features in the models. Another Ca x line pair (AA 1462, 
1504 A) is possibly present in absorption in NGC 246 and sug- 
gests a roughly solar Ca abundance. The only other object in 
which we discovered the Cax AA 1137, 11 59 A emission lines 
is the [WCE]-type central star NGC 2371. This corroborates the 
extraordinarily high effective temperature of this object. 
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